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We study the low- frequency, long- wavelength dynamics of liquid crystal elastomers, crosslinked 
in the smectic-j4 phase, in their smectic-yl, biaxial smectic and smectic-C phases. Two different 
yet related formulations are employed. One formulation describes the pure hydrodynamics and 
does not explicitly involve the Frank director, which relaxes to its local equilibrium value in a 
non-hydrodynamic time. The other formulation explicitly treats the director and applies beyond 
the hydrodynamic limit. We compare the low-frequency, long-wavelength dynamics of smectic-A 
elastomers to that of nematics and show that the two are closely related. For the biaxial smectic and 
the smectic-C phases, we calculate sound velocities and the mode structure in certain symmetry 
directions. For the smectic-C elastomers, in addition, we discuss in some detail their possible 
behavior in rheology experiments. 
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I. INTRODUCTION 

Liquid crystal elastomers [l| are soft amorphous solids 
that have the macroscopic symmetry properties of liq- 
uid crystals [2| . Usually, they consist of crosslinked side- 
chain or main-chain liquid crystal polymers. In their 
smectic phases, these elastomers possess a plane-like, 
lamellar modulation of density in one direction. In the 
smectic- A (SnxA) phase, the director describing the av- 
erage orientation of constituent mesogens is parallel to 
the normal to the smectic layers whereas in the smectic- 
C (SmC) phase, it has a component in the plane of the 
layers. 

Various unusual properties of smectic elastomers have 
been discovered to date experimentally and/or theoreti- 
cally, some of which are truly remarkable. For example, 
in smectic elastomers the crosslinking can suppress the 
Landau-Peierls-Instability [H, 3, HI that leads to a break- 
down of one-dimensional long-range order in other lay- 
ered systems such as stacked surfactant membranes (f| 0] 
or conventional smectic liquid crystals [8| . Depending on 
the chemical composition of the material 0, the stabi- 
lizing effect of the elastic network can be strong enough 
to allow for the persistence of one-dimensional long-range 
order over several microns fiol ] . SmA elastomers exhibit a 
non-Hookean elasticity when stretched along the normal 
of the smectic layers with a high Young's modulus Yt\ for 
strains below a threshold of about 3% and a consider ably 
smaller Yii for strains above that threshold [ill, [3 Il3j . 
SmC elastomers prepared by crosslinking in the SmA 
phase followed by a cooling into the SmC phase, like 
nematic elastomers 0, [HI, [l(| [I?], EH , are predicted to 
exhibit the phenomenon of soft elasticity 0, H3, [5l| . The 
phase transition from the SmA to the SmC phase is as- 
sociated with a spontaneous breaking of the rotational 
invariance in the smectic plane. The Goldstone mode of 
this symmetry breaking is manifests itself in the vanish- 



ing of certain shear moduli as well as the vanishing of the 
energy cost of certain nonlinear extensional strains. 

It is the interplay between elastic and liquid crystalline 
degrees of freedom that leads to these and other phenom- 
ena that are neither present in conventional rubber nor 
in conventional liquid crystals. The effects of this in- 
terplay are strongest in ideal, monodomain samples of 
liquid crystal elastomers, with a homogeneous orienta- 
tion of constituent mesogens. Generically, liquid crystal 
elastomers tend to segregate into many domains, each 
having its own local Frank director n specifying the di- 
rection of mesogenic order. In order to avoid such poly- 
domain samples, elaborate synthetization methods have 
been developed and applied to smectic elastomers for 
more than ten years now [22|, |23j, [2J| . Most fruitful, per- 
haps, have been multistage crosslinking techniques that 
allow for monodomain samples that neither need to be 
confined between director-aligning plates nor need to be 
kept in director-aligning external fields. These techniques 
have been used, for example, to synthesize elastomers 
with a permanent, macroscopically ordered SmA mon- 
odomain structure [25j as well as monodomain SmC elas- 
tomers that exhibit reversible phase transitions to the 
SmA and isotropic phases in response to temperature 
change [H, H3- 

Despite the successes in understanding the static prop- 
erties of smectic elastomers and the advances in their syn- 
thesis, dynamical experiments on these materials, such as 
rheology experiments of storage and loss moduli or Bril- 
louin scattering measurements of sound velocities, have 
not been reported to date. To interpret these kinds of ex- 
periments once they are done, it is desirable to have the- 
ories for the low-frequency, long-wavelength dynamics of 
smectic elastomers. Recently, we introduced such a the- 
ory for SmC elastomers [30( . Here, we present this work 
in more detail and we extend it to the SmA and biaxial 
smectic phases. To keep our discussion as simple as pos- 
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sible, we will exclusively consider elastomers, that have 
been crosslinked in the SmA phase, so that the smectic 
layers are locked to the crosslinked network [28|, H^] . To 
be specific, we will consider SmA elastomers crosslinked 
in the SmA phase and soft SmC elastomers that form 
spontaneously from these Srm4 elastomers when temper- 
ature falls below the SmA-to-SmC transition tempera- 
ture. In addition to these materials, we study soft biax- 
ial smectic elastomers that also could form spontaneously 
from a SmA elastomer, at least in principle, but that may 
be, as biaxial phases in general, hard to find in nature. 
We include these elastomers in our study mainly because 
their dynamics is similar to but simpler than that of soft 
SmC elastomers and, therefore, has a potential pedagog- 
ical value. 

Our theories fall into two categories. The first category 
is pure hydrodynamics which describes the dynamics of 
those degrees of freedom whose characteristic frequencies 
u> vanish as wavenumber q goes to zero. These theories 
focus on the very leading low-frequency, long-wavelength 
behavior and apply for frequencies u such that lot <C 1, 
where r is the longest non-hydrodynamic decay time in 
the system. As in conventional elastic media, our hydro- 
dynamic theories involve only the elastic displacement 
field u and not the director n, which relaxes to the lo- 
cal strain in a non-hydrodynamic time r n . The second 
category of theories explicitly includes n and applies for 
a larger range of frequencies than pure hydrodynamics. 
One motivation for setting up these theories is to bet- 
ter understand the role of n in the dynamics of smectic 
elastomers. The other motivation is that dynamical ex- 
periments, like rheology measurements, typically probe a 
wide range of frequencies that extends from the hydro- 
dynamic regime to frequencies well above it, and that, 
therefore, theories going beyond hydrodynamics could be 
valuable for interpreting these experiments. Our theo- 
ries with displacement and director assume that r n is 
the longest non-hydrodynamic relaxation time, r = t„, 
and that r„ is well separated from the next longest non- 
hydrodynamic time te, which we refer to as elastomer 
time. The specific origin of te is not essential to our 
work. For example, te might stem from Rouse-like dy- 
namics of the polymers constituting the elastomeric ma- 
trix. Provided the assumption r„ 3> te holds, our dy- 
namics with u and n applies for oj r„ <C 1 and lot„ > 1 . 

The outline of the remainder of this paper is as follows. 
Section [TT] briefly reviews the well established Poisson 
bracket formalism for coarse grained variables as well as 
the resulting equations of motion, when this formalism is 
applied to liquid crystal elastomers. Section UTT1 discusses 
the dynamics of SmA elastomers with strain and director. 
Sections IIVI and |V] treat, respectively, the hydrodynam- 
ics of biaxial smectic and SmC elastomers. Section IVT1 
presents our theory for SmC's with strain and director. 
Section [VIII contains concluding remarks. There are two 
appendixes that, respectively, compile results for elastic 
constants and provide some details of our calculations. 



II. DYNAMICAL EQUATIONS FOR LIQUID 
CRYSTAL ELASTOMERS 



In Ref. [3l| , where we studied the dynamics of nematic 
elastomers, we applied the Poisson bracket formalism for 
coarse-grained variables to derive hydrodynamical equa- 
tions for the displacement and momentum density as well 
as equations of motion that in addition involve the Frank 
director. The validity of these equations is, however, not 
limited to nematics. Our hydrodynamical equations ap- 
ply for conventional elastomers as well as for the usual 
(one component) liquid crystal elastomers in, e.g., their 
isotropic, nematic and smectic phases. Our dynamical 
equations involving n are valid for any liquid crystal elas- 
tomer with a defined director, like nematic and smec- 
tic elastomers. In this section we review the dynamical 
equations derived in Ref. [3l[ to provide important back- 
ground information, to stage known results that we will 
need as input as we move along and to establish notation. 



Poisson bracket formalism for coarse-grained 
variables 



Stochastic dynamical equations for coarse-grained 
fields [13, HH can be obtained by combining the Poisson- 
Bracket formalisms of Classical mechanics [34| and the 
Langevin 35] approach to stochastic dynamics. The 
Poisson-Bracket formalism ensures that the resulting 
equations of motion are, in the absence of dissipation, 
invariant under time reversal. The Langevin approach 
provides a description of dissipative processes and noise 
forces. Let ^(x, t), fi — 1, 2, ... be a set of coarse-grained 
fields that describe the long-wavelength, low-frequency 
dynamics of a system, i.e., they are hydrodynamic or 
quasihydrodynamic variables whose characteristic decay 
times in the long- wavelength limit are much larger than 
microscopic decay times. The dynamical equations for 
$ p (x,i) are first-order differential equations in time: 



Jd d x'J dt'{%(x,t),^>,t')} § ™ 



5H 



(x',f) 

(2.1) 



where TL is the a coarse-grained Hamiltonian that de- 
scribes the statistical mechanics of the fields "^(x, t). In 
Eq. (|2.1[) and in the remainder of this paper the Einstein 
summation convention is understood. The first term on 
the right hand side is a non-dissipative velocity, also 
known as the reactive term, that contains the Poisson 
bracket {$ M (x, i), (x', t')} of the coarse grained fields. 
The reactive term couples $ M to 5Ti/5<& v only if $ M and 
$„ have opposite signs under time reversal (when exter- 
nal magnetic fields are zero). The second term on the 
right hand side is a dissipative term. is the so-called 
dissipative tensor. It may depend on the fields $ M (x, t) 
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and it may contain —didj, where di = d/dxi. Its specifics equations of motion 
are determined by three principles: (i) r Mi „ couples $ M 
to STt/S^L, only if and <!>„ have the same sign under 
time reversal, (ii) By virtue of the Onsager principle [3^ | , 
it must be symmetric in the absence of external magnetic 
fields, (iii) T^^ has to be compatible with the symme- 
tries of the dynamic system. 

When Eq. (|2.1[) is augmented by a noise term C^, 
it represents a stochastic or Langevin equation that 
cound be used, for example, to set up a dynamic func- 
tional [13, [H, H^] to study the effects of nonlinearities 
and fluctuations via dynamical field theory. In this pa- 
per, we will ignore the effects of noise. 



1 

= ~9i, 
P 

9i = -7 h Vijki d 3 di u k 



(2.2a) 
(2.2b) 



B. Hydrodynamic equations 

Pure hydrodynamics describes the leading low- 
frequency, long- wavelength behavior of a system, i.e., it 
accounts exclusively for those degrees of freedom whose 
characteristic frequencies lu vanish in the limit of vanish- 
ing wavenumber. There are two general classes of hy- 
drodynamic variables: conserved variables and broken 
symmetry variables. A single-component liquid crystal 
elastomer has five conserved variables, viz. the energy 
density e, the mass density p, and the three components 
of the momentum density g. There are three broken- 
symmetry variables, namely the three components of the 
elastic displacement. Throughout this paper, we will use 
Lagrangian coordinates in which x labels a mass point 
in the unstretched (reference) material and R(x) labels 
the position of the mass point x in the stretched (tar- 
get) material, so that the broken-symmetry displacement 
variable is u(x) = R(x) — x. Since liquid crystal elas- 
tomers are, as any elastomers, permanently crosslinked, 
the mass density and the displacement are not indepen- 
dent. Rather, changes Sp in mass density are locked 
to changes in volume such that in the linearized limit 
Sp/p = — V • u. Thus, there are in total 7 independent 
hydrodynamic variables and 7 associated hydrodynamic 
modes. 

In the following, we will consider only isothermal pro- 
cesses so that the energy density, which is associated with 
heat diffusion, can be ignored. This leaves us with six hy- 
drodynamical variables, the components of the momen- 
tum density <7i(x) and the displacement u,(x). These 
are independent variables at the reference point x that 
satisfy the continuum generalizations of the usual rela- 
tions for the momentum and displacement of a parti- 
cle: Sgi(x)/Sgj(x') — 6ijS(x — x'), Sui(x)/6uj(x') = 
5ijS(x — x'), and Sgi(x)/Suj(x r ) = 0. These relations 
yield a non- vanishing Poisson bracket between Ui (x) and 
.gj(x'): {ui(x),gj(x')} = S lj S(x - x'). The g - g and 
u u Poisson brackets are zero. Using these results 
and the fact that the coarse-grained kinetic energy reads 
"Hkin — J d 3 x g 2 /(2/o), one obtains the coarse-grained 



Here an in the following we use the convention that in- 
dices from the middle of the alphabet, {i,j, k, I}, run 
from 1 to 3 (corresponding to the x, y, and z directions). 
H = Jd 3 xf is the elastic energy describing the elas- 
tomer with / the corresponding elastic energy density. 
The specifics of / depend on the particular elastomeric 
phase under consideration and will be discussed as we go 
along. The absence of any dissipative term proportional 
to —SH/Sui in Eq. (|2.2aj) . which would describe perme- 
ation, reflects the tethered or crosslinked character of the 
elastomer. Pure couplings of iji to ilk are forbidden be- 
cause Eqs. (|2.2I) have to obey Galileian invariance, and 
thus the components of the dissipative tensor coupling 
cji to iik are of the form —rjijki djdi, with rjijki being the 
viscosity tensor. In our discussion of the hydrodynamics 
of the biaxial smectic and the SmC phase, to be pre- 
sented in sections IIVI and [V] we will restrict ourselves for 
simplicity to frequencies uj that are much less than the 
inverse characteristic time r^ 1 associated with the vis- 
cosities. In this limit, the rj^ki can be viewed as local in 
time, or equivalently, their temporal Fourier transforms 
can be considered as being constant. 

To discuss the dynamics of the individual smec- 
tic phases in detail, it will be convenient to combine 
Eqs. (|2.2ap and (|2.2b[) into a set of second order differ- 
ential equation for the displacement components Ui only. 
Switching from time to frequency space via Fourier trans- 
formation, this set of equations can be expressed as 



pu> 2 Ui 



(2.3) 



where cr^ (tu) are the components of the stress tensor a. 
The specifics of the stress tensor, which depend on / 
and the viscosity tensor, will be stated for the individual 
phases as we proceed. 



C. Dynamic equations with displacement and 
director 

When the Frank director is included as a dynamical 
variable, there is an additional non-vanishing Poisson 
bracket, viz. {rii(x), gj(x')} = -X ljk dk5(x- x') [I3,[iH. 
Then, Eq. (|2.ip leads to the equations of motion 



■ fJ, 



5H 

5m ' 



SH 5H 
Ajik Cfe- h i>ijki OjOiUk 

OUj OUi 



(2.4a) 
(2.4b) 
(2.4c) 
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The properties of the tensor A^fc are dictated by three 
constraints: First, the magnitude of the director has to 
be conserved, i.e., n ■ n = implying riiXijk = 0; second, 
the equations of motion must be invariant under n — > 
— n implying A^fe must change sign with n; and third, 
under rigid uniform rotations, the director has to obey 
n = |(V x u) x n. The constrains imply that Xijk is of 
the form 

hjh = ^ ( 5 ij n k + S 7k n j) - \ ( S fj n k ~ Sfk nj) , (2.5) 

where Sfj = 5ij — riirij is the projector on the subspace 
perpendicular to n. There are no constraints on the value 
of A. The second term on the right hand side of Eq. (|2.4a[) 
is a dissipative term that describes diffusive relaxation 
of the director. Its coefficient T has the dimensions of 
an inverse viscosity. As before, there are no dissipative 
contributions to Eq. (12.4b[) because liquid crystal elas- 
tomers are tethered. There is no dissipative contribution 
to Eq. (|2.4a[) involving 5H / Sui because the terms on the 
right side of the equation must be odd in nt owing to 
n > — n symmetry. Any such terms must be of the 
form Aij STC/Suj, and it is impossible to construct a ten- 
sor Aij that satisfying the symmetries of the system and 
making STi/Suj odd under n — > — n. A similar ar- 
gument rules out a dissipative contribution to Eq. ()2.4b[) 
that involves SH/Snj. As before, Galileian invariance de- 
mands that the dissipative coupling of gi to ilk are of the 
form —Vijki djdi. Here we use the notation Vijki instead 
of rjijki for the viscosity tensor, so that we can cleanly 
keep track of differences between purely hydrodynamical 
theories and theories with strain and director. In gen- 
eral, we will assume that T and the i/yjw are, like the 
rjijki, independent of frequency. In Sec. IVI C[ where we 
discuss the behavior of SmC elastomers in rheology ex- 
periments, however, we will refrain from this simplifying 
assumption and we will replace the constants T and the 
Vijki by phenomenological functions T(uj) and Vijki{io) of 
frequency. 

III. SMECTIC-A ELASTOMERS - DYNAMICS 
WITH STRAIN AND DIRECTOR 

Macroscopically, SmA elastomers have uniaxial sym- 
metry. When they are crosslinked, as we assume, in the 
Sim4 phase, the uniaxial symmetry is not spontaneous 
but rather permanently imprinted and the material can- 
not possess soft elasticity. In the small deformation limit, 
these elastomers are macroscopically simply uniaxial rub- 
bers (albeit with a very large value of the modulus C\ 
for extension or compression along the layer normal) . As 
such, their hydrodynamics is identical to that of uniaxial 
solids 42| Beyond the hydrodynamic regime, however, 
one must expect an influence of the director on the dy- 
namics of SmA elastomers in which the genuine differ- 
ence between SmA and conventional uniaxial elastomers 
shows up. Here we apply our dynamical equations with 



displacement and director to the SmA phase in order to 
study these effects. 



A. Elastic energy 

As alluded to above, SmA elastomers are macroscop- 
ically simply uniaxial rubbers, at least when strains are 
small. Therefore, the harmonic stretching elastic energy 
density / u of a SmA elastomer has exactly the same form 
as that of a uniaxial solid [43| : 

/ u = \ Cx u 2 zz + C 2 u zz u u + \ C 3 u% 

+ C 4 ul b + C 5 u 2 az . (3.1) 

We omit higher order terms in the strain because they are 
inconsequential for our linearized dynamical equations. 
We have chosen the coordinate system so that the z-axis 
lies in the uniaxial direction. Here and in the following, 
indices from the beginning of the alphabet, {a, b}, take 
on the values 1 and 2. The Uij are components of the 
Cauchy-Saint-Venant 0, |45| strain tensor u, 

uy(x) = £(A£Ajy -Sij) (3.2a) 
= \ {d iUj + djUi) + • • • , (3.2b) 

where Ay = dRi/dxj are components of the Cauchy de- 
formation tensor A. In Eq. (|3.2bp we have dropped the 
nonlinear part of the strain tensor, because, for our pur- 
poses, it is sufficient to work to harmonic order. 

For setting up a dynamical theory for Sim4 elastomers 
that goes beyond the hydrodynamic limit, we need a 
more detailed level of description that explicitly involves 
the Frank director n. When n is included, the overall 
clastic energy density / will have the following contribu- 
tions: 

/ = /u + /n-tilt + /n-Frank + /coupl : (3-3) 

with /u as given in Eq. (|3.ip . / n -tnt is a tilt energy den- 
sity that accounts for the energy cost when the director, 
which in an equilibrium SmA elastomer is aligned along 
the normal N of the smectic layers, tilts away from N. 
fa- Frank is the density of the usual Frank energy that 
describes deviations from the homogeneous equilibrium 
orientation of the mesogenic component. / COU pi finally is 
the energy density of coupling between strain and direc- 
tor distortions. 

Coupling the strain and the director brings about an 
intricate conceptual problem: whereas the liquid crys- 
talline fields N and n transform as (rank 1) tensors in 
target space and are scalars with respect to rotations in 
reference space, the strain tensor u is a (rank 2) tensor in 
reference space, and it is a scalar in target space. Thus, 
in order to construct meaningful combinations of u and 
the liquid crystalline fields, we have to de able to repre- 
sent these quantities in either space. The matrix polar 
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decomposition theorem [46[ provides a route to this rep- 
resentation. It allows us convert (or rotate) any reference 
space vector b to a target space vector b via 

b = Qb (3.4) 

and a target-space vector to a reference space vector via 

b = Q T h, (3.5) 

with the ration matrix O given by 

2 = 4(4 T 4r 1/2 =l + R A + --- > ( 3 - 6 ) 

where 1 is the unit matrix and where r\ ^ = \(r\ — r] T ) is 
the antisymmetric part of the displacement gradient ten- 
sor r\ defined by rjij = djUi. In particular, we can rotate 

the target space director n = (c, n z ) , n z = yl — cf, 
where c is the so-called c-director, to a reference space 
vector n = (c, n z ) , fi z — -^/l — c|, and rotate the unit 
vector e = (0, 0, 1) specifying the direction of the uniaxial 
axis in the reference space to a target-space unit vector. 

With the use of the polar decomposition technique, it is 
not difficult to formulate the tilt and the coupling energy 
in reference space variables. As mentioned above, in the 
equilibrium SmA phase, the director prefers to be paral- 
lel to both the layer normal N and the anisotropy axis 
e (47[, which are parallel to each other. There are energy 
costs proportional to (N • n) 2 and (e • n) 2 associated with 
deviations from this equilibrium. These combine to yield 
a contribution to the tilt energy proportional to c 2 and 
higher order terms that are unimportant here. Therefore, 
the leading contribution to the tilt energy density is 

/„-ti 1 t = ir C . 2 = irQ 2 + ... , (3.7) 

where Q a — c a — r\Aaz- At harmonic order, there is only 
one contribution to the coupling energy that is compati- 
ble with the uniaxial symmetry, viz. 

/coupi = A4 c a u az = A4 Q a u az + ■ ■ ■ ■ (3.8) 

(the symbols Ai, A2 and A3 are usually reserved for 
third-order couplings that, however, do not matter here). 
When /n-tiit and / CO upi are combined, the result is iden- 
tical in form to the harmonic coupling energy density 

/ u ,n = \ D x Ql + D 2 u za Q a , (3.9) 

appearing in the elastic energy density 

fnem = fu ~^ /n— Frank + /u,n , (3.10) 

for a nematic elastomer, when the latter is described in 
terms of strains and the Frank director [48| . The stretch- 
ing energy density / u and the Frank energy density 
/n-Frank in Eq. (|3.10p are the same as in Eq. (|3.3p . albeit 
with the elastic constant having, in general, different val- 
ues. Moreover, r and A4 will be different from D\ and D 2 . 



In nematic elastomers, the values of D\ and D 2 are such 
that the shear modulus C5 is renormalized by director 
relaxation to a value C^_ ncm = C5 — D\jD\ that is zero 
for an ideally soft nematic elastomer or nearly so for a 
non-ideal semi-soft nematic elastomer. The correspond- 
ing renormalized shear modulus C^_ SmA = C5 — A|/r 
of a Sm^4 elastomer, to the contrary, will always be 
significantly larger than zero, when the elastomer was 
crosslinked, as we assume, in the SmA phase. 

B. Low- frequency, long- wavelength dynamics 

Inserting the just discussed elastic energy of SmA elas- 
tomers into the equations of motion with displacement 
and director, Eq. (|2.4p . we obtain dynamical equations 
for Srm4's that are identical in form to those of nematic 
elastomers. The only differences reside in the magnitude 
of the elastic constants, viscosities and nonhydrodynamic 
relaxation times. In addition, C§_ nom can vanish in soft 
nematics but C^_ SmA is always significantly larger than 
zero in SmA elastomers crosslinked in the SmA phase. 
Thus, solving the SmA equations leads to sound veloc- 
ities, modes etc. that have exactly the same from as in 
nematics. To avoid undue repetition, we will refrain from 
reviewing these results here in detail; rather we refer the 
reader directly to Ref. [3lj . 

IV. BIAXIAL SMECTIC ELASTOMERS - 
HYDRODYNAMICS 

When the shear modulus C4 in the elastic energy den- 
sity (|3.ip becomes negative, as it will in response of 
biaxial ordering of the constituent mesogens of a SmA 
elastomer, the system undergoes a phase transition to 
a broken-symmetry state with D 2 h (orthorhombic) sym- 
metry [l^ . [20| . This mechanism, at least in principle, 
could pro duce a soft or semi-soft biaxial smectic elas- 
tomer |l9l . I20I l2lj . In this section we study the dynam- 
ics of these elastomers. In the framework of Lagrange 
elasticity theory, a model for their elastic energy den- 
sity exists thus far in a strain-only formulation but not 
in a formulation with strain and director. Thus, we will 
restrict ourselves here to pure hydrodynamics. 

A. Elastic energy, viscosity tensor and 
hydrodynamical equations 

A biaxial smectic elastomer is macroscopically an elas- 
tic body with D 2 h symmetry. Therefore, its stretching 
energy is that of orthorhombic systems [431 ]. albeit with 
an important difference: the modulus for shears in the 
plane of the smectic layers (with our choice of coordinates 
the ly-plane), C^ yxy , vanishes for an ideally soft biaxial 
smectic elastomer or is small for a semi-soft one. For de- 
tails on the derivation of the stretching energy density of 
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soft biaxial smectic elastomers, we refer to Refs. [19ll20l|. 
If C^ y vanishes, shears u xy in the xy-plane cost no elas- 
tic energy and, therefore, cause no restoring forces. Thus, 
bending terms proportional to (d 2 u x ) 2 and (d 2 u y ) 2 have 
to be added to the harmonic elastic energy to ensure me- 
chanical stability. This leads to a overall elastic energy 
density of the form 

/ = 2 C ZZZZ U 2 ZZ + i C XZXZ U xz + i CyzyZ Uy Z 

~T" C ZZXX U ZZ U XX "t" C ZZ yy U zz Uyy ~\~ ^ C XXXX V> xx 

2 Cyyyy U yy + C xxyy U xx U yy + C xyxy U xy 

+ l Bl {dlu x ) 2 + \B 2 {dlu y )\ (4.1) 

with bending moduli B\ and B 2 . 

Like the elastic constant tensor ki > the viscosity 
tensor rjijki of an orthorhombic system has 9 indepen- 
dent parameters. It can be parameterized so that the 
entropy production density Ts takes on the same form 
as the stretching energy density, with Cijki replaced by 
rjijki and u%$ replaced by the linearized form of Uij, i.e. 
Uij = diUj + djiii, 

Ts = J Vzzzz U zz + J Vxzxz u X z + 5 Vyzyz u yz 

1 • 2 

~t~ IJzzxx ^zz^xx Vzzyy ^zz^yy 2 Vxxxx ^xx 

2 Vvyyy ^yy Vxxyy u xx u yy + ^ r/ xyxy u xy . (4.2) 

Equation (|4.2p contains all contributions to Ts in the hy- 
drodynamic limit. In writing it, we have assumed that 
all non-hydrodynamic degrees of freedom, like the direc- 
tor, have relaxed to their local equilibrium values in the 
presence of My . r] xyxy is an effective shear viscosity that 
is, like the corresponding shear modulus C xyxy , renor- 
malized by the relaxation of the director. 

Now we possess all the ingredients that are required 
to write down the hydrodynamical equations for biaxial 
smectics. Inserting the elastic energy (|4.1[) into Eq. (|2.2p 
and using the form of the viscosity tensor implied in 
Eq. (|4.2[) , we obtain the equations of motion in the form 
of Eq. (|2.3[) with the components of the stress tensor g_ 



jiven by 








&xx{oj) = 


Cxxxxil-tj) ^xx ~t" 


Cxxyyi^*) V*yy 


+" C XXZ z{u) U zz , 

(4.3a) 


Vyyi^) = 


C xxyy {iO^) u xx -f- 


Cyyyyi 1 ^) u yy ' 


h Cyy ZZ {uj) Uzz , 

(4.3b) 


<Jzz{u) = 


C xxzz (cj) u xx + 


Cyyzz^) u yy ^ 


■ C zzzz (uj) u zz , 
(4.3c) 


<?xy{u) = 


2 ^xyxyi 10 ) u xy 


- Bid^u x , 


(4.3d) 


&yx(u) = 


\ C X y X y(u) U X y 


- B 2 dlu v , 


(4.3e) 


O-xz(u) = 


& zx{u) — \ C xzxz {uS) u xz , 


(4.3f) 


Vyz{u) = 


0Z V {U) = \ Cy Z yz{u)Uyz ■ 


(4.3g) 



Here, Cijki{uj) stands for 



C^ki (w) = Cijki - iw mjki , (4.4) 



and the renormalized frequency dependent modulus 
C xyxy (uj) is defined as 

^xyxyi.^*) ^ xyxy Vxyxy ■ (4-5) 

When the renormalized shear modulus C xyxy is nonzero, 
in which case the bending terms featuring B\ and B 2 are 
unimportant for the hydrodynamic behavior, the stress 
tensor (|4.3|) is identical to the hydrodynamic stress tensor 
for an orthorhombic solid and therefore the equations of 
motion ()2.3|) are dentical to the hydrodynamic equations 
for an orthorhombic solid in this case. 



B. Sound velocities 

To assess the mode structure of biaxial smectic elas- 
tomers, we start with an analysis of propagating sound 
modes in the dissipationless limit, i.e., in the limit where 
viscosities rjij^i are zero. The sound modes have frequen- 
cies 

w(q) = C(#,<p)q, (4.6) 

where q = |q| and where and ip are the azimuthal 
and polar angles of the wavevector q in spherical coor- 
dinates. In calculating the sound velocities, one can ne- 
glect the bending terms in the stress tensor (|4.3p . even 
though these are nondissipative. As we will see further 
below, the bending terms give rise to modes with fre- 
quencies to ~ q 2 along the symmetry directions, where 
sound velocities vanish, and these modes mix with dissi- 
pative ones to become overdamped diffusive modes with 
to ~ —iq 2 . Therefore, the bending terms do not con- 
tribute to the sound velocities of the propagating modes. 
For simplicity, we will focus here on the ideally soft case 
with C xyxy = 0. We will return to the more general case, 
where C xyxy can be nonzero, in Sec. lIVCl 

Setting the viscosities, the bending moduli K± and K 2 
and the shear modulus C xyxy to zero and switching from 
the reference space coordinate x to the wavevector q via 
Fourier transformation, the hydrodynamical equations 
simplify to 

p io il x — \Cxxxx q x ~l~ 4 C XZX z q z ~\ V'x C XX y y q x qy V> y 

~l~ \_C XXZZ ~~\~ t C xzxz ^q x q z u z , (4.7a) 

P UJ Uy — Cxxyy QxQy ^x ~^ \_Cyyyy Qy H~ 4 Cy Z y Z Cf z ~^ Uy 

+ \Cyyzz + \ Cyzyz] q y q z u z , (4-7b) 

p UJ U z — [Cxxzz ~t~ ^ C X zxz\QxQz ^X 
H~ \Cyyzz 4 Cy Z y Z ~^ QyQ z Uy 

+ [4 C XZXZ q x + J Cy Z y Z Qy + C ZZ ZZ Q Z ] U Z . 

(4.7c) 

Inserting the frequency (|4.6[) into Eq. (|4.7p and express- 
ing the components of q in spherical coordinates, one 
can solve for the sound velocities C{d,<p). The result- 
ing sound velocities are plotted schematically in Fig. [T] 
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FIG. 1: Schematic plots of sound velocities, (a), (b) and (c): 
spherical plots of respectively mode pairs (i), (ii) and (iii). (d) 
and (e): polar plots in the xy- and zz-planes, respectively, of 
all 3 sound-mode pairs. 



There are 3 pairs of sound modes. One of these pairs (i) 
is associated with soft shears in the xy-plane. Its veloc- 
ity vanishes for q along the x and y directions, so that 
when viewed in the xy-plane it has a clover-leave shape. 
The remaining 2 pairs are associated with non-soft de- 
formations. In the incompressible limit, these pairs be- 
come purely transverse (ii) and longitudinal (iii), respec- 
tively. In the xy and irz-plancs, their velocities are non- 
vanishing in all directions. 

Figure Q] indicates that biaxial smectic elastomers, if 
they can be produced, have the potential for technolog- 
ical application. Since the mode pair (iii) is purely lon- 
gitudinal in the incompressible limit, its velocity will be 
much greater than that of the remaining two pairs. The 
sound velocities of the remaining pairs, which have differ- 
ent transversal polarizations, are so that the velocity of 
pair (i) vanishes in directions where the velocity of pair 
(ii) remains finite. Therefore, a biaxial smectic elastomer 
could in principle be used to separate the sound mode 
pair (i) from the sound mode pair (ii), or in other words, 



such an elastomer could in principle be utilized in mak- 
ing acoustic polarizers. This kind of application has been 
proposed a while ago for soft nematic elastomers [HH • In 
our discussion of the hydrodynamics of soft SmC elas- 
tomers, Sec. El we will see that SmC's also allow splitting 
of sound waves. It seems, as if the potential for acous- 
tic polarization is a generic feature of soft liquid crystal 
elastomers. 



C. Mode structure in the incompressible limit 

Now we turn to the mode structure, including dissi- 
pation, in the incompressible limit. Even in this limit, 
the equations of motion remain fairly complicated and 
their solutions have the potential for unpleasant alge- 
braic complexity. Therefore, we will restrict ourselves in 
the following to the two simplified cases, where either q y 
or q z vanishes, as respectively in Fig. Q] (d) and (e). De- 
spite these simplifications, the resulting equations of mo- 
tion are still to complicated to be solved in closed form. 
However, since we are interested in the long-wavelength 
behavior, we can perturbatively determine solutions in 
the from of a power series in the wavevector. These so- 
lutions will turn out to be of the following two types: (i) 
propagating modes with frequencies 

oj p = ±Cq-iD p q 2 (4.8) 

with sound velocities C and diffusion constants [4i| D p , 
and (ii) diffusive modes with frequencies 

uj d = -iD d q 2 ± ^-{D d q*y- + Bq± (4.9) 

with diffusion constants D d and bending terms B. For 
B /D 2 <C 1 the diffusive modes split up into slow and fast 
modes 

uj s = -iB/{2D d )q 2 , (4.10a) 
uj f = -i2D d q 2 . (4.10b) 

Specifics of the sound velocities, the diffusion constants 
and the bending terms are given in the following. 
For q y = 0, the equations of motion reduce to 

pu 2 u x = [C xxxx (lj) q 2 x + \ C xzxz (o>) q 2 ] u x 

+ [C xxzz (uj) + jC xzxz (uj)]q x q z u z , (4.11a) 

pu 2 u y = [\C* yx y(u)ql + \C yzyz {uj) q 2 + B 2 qt]u y , 

(4.11b) 

puj 2 u z = [C xxzz (lj) + jC xzxz (uj)]q x q z u x 

+ [§ C xzxz {uj)q 2 x + C zzzz (uj)q 2 z ]u z . (4.11c) 

One of the simplifications that we can enjoy for q y = 
is that the equation of motion for u y decouples from the 
equations of motion for u x and u z , which are coupled to 
each other. The u y equation produces a pair of transverse 
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propagating modes polarized in the y-direction, u 
with a sound velocity 



with q z replaced by q y as well as C xyxy and rj^ replaced 
by C xzxz and t] xzxz , respectively: 



Cy — 



^xyxy Qx + C yzyz q z 



4p 



and a diffusion constant 



D 



Vxyxy Qx + Vyzyz Q z 



p,y 



(4.12) 



(4.13) 



where qi stands for qi/q. For the ideally soft elastomer, 
Cxyxy = 0' the sound velocity vanishes for q z = 0, and the 
above modes become diffusive with a diffusion constant 



p, 'Ixyxy ~2 



and a bending related constant 

B y = Q4pB 2 . 



(4.14) 



(4.15) 



In the limit 64pB 2 (Vxyxy) 2 > these diffusive modes for 
q z — split up into a slow and a fast diffusive mode with 
frequencies 



45 



2 „2 



Vxyxy 



Q x > 



Vxyxy 9 
-l qt ■ 



4p 



(4.16a) 
(4.16b) 



The equations of motion for u x and u z can be solved 
by decomposing (u x , u z ) into a longitudinal part ui along 
q and a transversal part ut- In the incompressible limit 
ui vanishes. The equation of motion for ut produces a 
pair of propagating transversal modes with polarization 
u || ex, where = (q z , 0, — q x ). This mode pair has the 
sound velocity 



C'i 



'4(0, 



C z 



2c xxzz )ej 2 z + c xzx M-f z y 



(4.17) 



C z 



D P,z = 



/ C xzxz q x + C yzyz q y 

V ~p 

TJxzxz (j x + Vyzyz <}y 



8p 



(4.19a) 
(4.19b) 



Since both elastic constants contributing to the sound 
C z are nonzero, even for an ideally soft biaxial smectic 
elastomer, this sound velocity is nonzero in any direction 
in the xy-plane. 

To solve the coupled equations for u x and u y , we pro- 
ceed as above, and we decompose {u x ,u z ) into a longi- 
tudinal part ui and a transversal part ut, which now is 
along e-r = (q y , 0, —q x )- In the incompressible limit, ui is 
suppressed, and of the two coupled mode pairs only the 
transverse pair polarized along survives. The sound 
velocity and the diffusion constant of this pair are, re- 
spectively, identical to Eqs. (|4.17p and (|4.18p . however 
with q z replaced by q yi with C ZZZZl C xxzz and C xzxz , re- 
spectively, replaced by C yyyy , C xxyy and C R yxy , and with 
corresponding replacements for the viscosities: 



Cn 



\4-{C XXXX + Cyyyy 2C 'xXyy)(}xQy + Cx\xy(Qx Cj y )' 



4p 



Vyyyy %Vxxyy ~2 ~2 , Vxyxy / ~2 _ ~2\2 
Hx^y ' o n ^"x Vy) 



2p 



(4.20a) 

~fy) 2 - 

(4.20b) 



Since now C xyxy appears in the sound velocity, the sound 
velocity vanishes for the ideally soft elastomer if q x = 
or q y = 0. In either case, the modes become diffu- 
sive. For q y = 0, the frequencies of these diffusive 
modes are identical to the frequencies of the diffusive 
modes discussed further above. Therefore, in the limit 
64 pK 2 <C (il xyxy ) 2 one again has a splitting into slow 
and fast diffusive modes with frequencies as stated in 
Eq (pTl"6]) . For q x = 0, we find 



and a diffusion constant 

^ Vxxxx IJzzzz 2tJ xxzz o ~2 

D p,t = ~ q x q z 

L p 



Ixzxz , x 2 -2\2 

Qz) ■ 



(4.18) 



We complete our discussion of the hydrodynamics of 
biaxial smectic elastomers by turning to the case that q 
lies in the xy-planc. For q z = 0, the equation of motion 
for u z decouples form those for u x and u y , which remain 
coupled to one another. The u^-equation yields a pair 
of transverse propagating modes polarized along e z with 
sound velocities and diffusion constants, which are of the 
same form as those given in Eqs. (|4.12p and (|4.13|) , albeit 



D d , T = 



Vxyxy _> 



i/0 



Qz, 



B T = 64p Bi 



(4.21a) 
(4.21b) 



In this case, a splitting up into slow and fast diffusive 
modes occurs for 64 pB\ <C (r] xyxy ) 2 . The frequencies of 
these slow and fast diffusive read 



Uf,T 



■ jgl 2 

1 n R q y 

Ixyxy 

■ Vxyxy o 

a Qy 
4p " 



(4.22a) 
(4.22b) 



9 



V. SMECTIC-C ELASTOMERS - 
HYDRODYNAMICS 

In this section we investigate the dynamics of SmC 
elastomers in the hydrodynamic limit. First we will re- 
view their elastic energy density in a model with strain 
and director [23]. This model will become particu- 
larly important in Sec. lVIl Here, we need a model elastic 
energy density in terms of strain only, which we derive 
from the more complete model by integrating out the di- 
rector. Then we will come to the actual hydrodynamics. 
We will write down the complete hydrodynamical equa- 
tions and we will extract the sound velocities and the 
mode structure. 



A. Elastic energy 

When formulated in terms of strains and the Frank 
director [2(| , the elastic energy density / of a SmC 
elastomer crosslinked in the SmA phase can be divided 
into three parts 



/ — /u + /u,n + /n— Frank ■ 



(5.1) 



where / u is the stretching energy density depending only 
on /u,n describes the coupling of the Frank director 
n to strain variable, and / n -Frank is the density of the 
Frank energy. In the following we choose the coordinate 
system so that the z-axis is parallel to the director of 
the initial SmA phase and the ir-axis is parallel to the 
direction of tilt in the resulting cqulibrium SmC phase. 
Stated more precisely, we choose our coordinates so that 
the equilibrium reference space director n° = (c°,n°) 
characterizing the undeformed SmC phase is of the form 
n° = (5,0,Vl - S 2 ), with S = £° x being the order pa- 
rameter of the transition. It should be emphasized that, 
in general, one has to distinguish carefully between the 
reference-space equilibrium director n° and the physical 
equilibrium director n°. Note from Eqs. (|3.4[) and (|3.6[) . 
however, that to leading order in the displacement gra- 
dients rjtj, n° and n° coincide, 



(S,0,y/l-S 2 ). 



(5.2) 



For our linearized dynamical theories, only that leading 
order matters and, therefore, for our current purposes, 
we can refrain from distinguishing n° and n°. 

With our choice of coordinates, / u can be written in 
the same form as the elastic energy density of conven- 
tional monoclinic solids [43| . 

/u = 2 Cxyxy U xy + C X yzy U xy U zy + ^ C zyzy U 2 y 



i 2 ^zzzz ^zz 1 2 ^xxxx "- X x 1 2 WW yy 
+ 2 Cxzxz U xz + C ZZX x U zz U X x + c z zyy UzzUyy 
~b Cxxyy Uxx^yy -f - Cxxxz ^xx^xz ~\~ Cyyxz ^yy^xz 
~\~ Czzxz Uzzl^xz ) (5.3) 



but with constraints relating the three elastic constants 
in the first row. These latter constants can be expressed 
in terms of an overall elastic constant C and an angle 
a which depends on the order parameter S, as C xy xy — 



Ccos 2 e, c x 



C cos sin 0, and C z 



CW0. 



The coupling energy density / Ujn can be stated as 



A, 



\ A [5c y + a u xy + (3 u yx ] 



(5.4) 



where A is a coupling constant, where a and f3 are di- 
mensionless parameters, and where 5c y is the deviation 
of c y from its equilibrium value c y = 0. In the following 
it will be more useful to work with the physical direc- 
tor than with its reference space counterpart. Switching 
from <5n = n — n° to Sn = n — n° via the transformation 
we get 



/ u ,n = 5 A [Q y + a u xy + (3 u yx f 
with the variable Q y defined by 



Qy 



5n,. 



S TjAyx - Vl - S' 2 



VAyz : 



(5.5) 



(5.6) 



and where higher order terms in rjij have been discarded. 

The remaining contribution to the total elastic energy 
density that we need to discuss here is the Frank energy 
density. When expanded to harmonic order about the 
equilibrium director (|5.2[) . it becomes 



/n— Frank 

— 2 ^xxxx i&xTix) ~t~ 2 -^yyyy ^y^y) o ^yxyx {.^x^iy) 

~\~ 2 -K-xyxy ij3y^x) o Kxzxz i^z^x) ~l~ o ^-yzyz (dzfly) 
~\~ K xxy y d x Tl x 9 y Tl y -(- K X yyx ^x^y &yTl x 
~\~ K xxxz d x Tlx d z 7l x ~\~ Kyyyx *3yTl y (3 y Tl x 

-(- K yxyz O x Tl y &z^y ~l~ -Kxyyz & y Tl x d z Tl y , (5.T) 

where we have replaced <5n by n for notational simplicity, 
with the understanding, that it has only two components 
n x and n y . We will stick to this abbreviated notation 
for the remainder of this paper. The elastic constants 
Kxxxx and so on are combinations of the the original 
Frank constants, describing respectively splay, twist, and 
bend distortions of the director and the order parameter 
S. For specifics, see appendix IA1I 

As discussed above, the director is not a genuine hy- 
drodynamic variable. Therefore, it should be integrated 
out of the elastic energy as long as we focus on pure hy- 
drodynamics. Minimizing / over n x and n y , we find that 
these quantities relax in the presence of strain to 



0. 



(5.8a) 



v = S rjAyx + \A - S 2 r)Ayz - a u xy - (3 1 



1 

2A 



Ky X yx X -\- Ky yyy dy -\- Kyxyz & z ~t~ ^Kyxyx (3x^3 z 



X I S TjAyx + - S 2 VAyz ~ au xy ~ (3 U yz } 



(5.8b) 
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where terms of fourth and higher orders in derivatives 
have been discarded. Inserting these results into /, we 
obtain an effective elastic energy density of the form 

/ = /» + /bend , (5.9) 

where / u remains as given in Eq. (|5.3p . and where /bend 
is a bending energy density given by 

/bend =\B 1 (dlUyf + \B 2 (d 2 z U y ) 2 + \ B 3 (d x d z U y ) 2 

+ B 4 u y d z dlu v + B 5 u y d x d z u y + \ B e (d 2 u x ) 2 
+ \ B 7 (d 2 y u z ) 2 + B 8 (d 2 u x )(d 2 u z ) . (5.10) 

Specifics of the bending moduli B\, B 2 and so on are 
given in appendix I A2I 

The elastic energy density (15 .9|) implies that SmC elas- 
tomers exhibit static soft elasticity. Suppose deforma- 
tions arc not too large so that the bending contribu- 
tion to Eq. (|5.9|) may be neglected. Then, due to the 
above relations among the elastic constants C xyxy , C xyzy 
and C zyzy , deformations characterized by Fourier trans- 
formed displacements u || e v and wavevectors q | e 2 = 
(— sin6>, 0, cos#) or, alternatively, u || e 2 and q || e y cost 
no elastic energy and hence cause no restoring forces. In 
the following, we will occasionally switch to a coordinate 
system that is rotated relative to the original coordinate 
system though an angle 9 about the y-axis so that axis of 
the new system are e x i — (cos 9, 0, sin#), e y > — e y , and 
e 2 ' = e 2 . One one hand, the effects of soft elasticity are 
more evident in the rotated system. On the other hand, 
the directions e x i and e z > will play a certain role in the 
mode structure of SmC elastomers, as we will see below. 
In the rotated system, the stretching energy density takes 
on the form 

/u 2 ^x'y'x'y' (^x'y') ~t~ 2 Cz'z'z'z' i^z'z 1 *) 
~t~ 2 ^x' z' x' z' i^x' z'*) "t" C z ' z ' x ' x ' 1L Z ' z'^x'x' 

~l~ C Z ' Z ' y' y> U z ' z l U y > y ' -)- ^ Cx'x'x'x' {^X'X'^} 

2 ^y'y'y'y' (.^y'y') dx'x'y'y' ^x' x'^y' y' 
~l~ C x ' x ' x ' z ' \L X > x 'U x > z i -\- C y ' y' x' z' Uy'y'U x ' z ' 
+ Cz'z'x'z' u z'z' u x'z' j (5-11) 

where C x > y > x ' y ' = C, C y > y > y ' y i = C yy yy and where the 
remaining new elastic constants are non-vanishing con- 
glomerates of the elastic constants defined via Eq. (|5.3p 
and sines and cosines of 9. Note that, in this coordinate 
system, the elastic constants C X ' y > y > z > and C y 'z' y 'z' are 
zero, and that therefore the elastic energy does not de- 
pend at all on u y i z i , impyling there is no energy cost for 
shears in the z'y'-plane. If the elastomer is crosslinked 
in the SmC phase, these moduli become nonzero. If they 
remain small, the elastomer will be semi-soft. 

B. Hydrodynamic equations 

Knowing the effective elastic energy density of SmC 
elastomers in terms of the displacement variable only, 



Eq. (|5.9p in conjunction with Eqs. (|5.3[) and (|5.10[) . we 
are almost in the position to write down hydrodynamic 
equations for SmC elastomers. The only ingredient that 
is still missing is the viscosity tensor. We will parametrize 
this tensor so that the density of the entropy production 
from viscous stresses takes on a form similar to that of 
Eq. (JO}: 

Ts = ^ Vxyxy ^xy Vxyzy ^xy^zy + 2 ^zyzy ^zy 

H~ 2 Vzzzz U zz + 2 Tjxxxx U xx + 2 Vyyyy ^yy 2 T l xzxz ^xz 

T]zzxx ^zz^xx \ Tjzzyy ^zz^yy H~ ^xxyy ^xx^yy 
~\~ Vxxxz ^xx^xz 1 TJyyxz UyyU X z 1 ^jzzxz ^zz^xz ■ 

(5.12) 

As in Eq. (|4.2[) , we assume here that all non- 
hydrodynamic degrees of freedom including the director 
have relaxed to their local equilibrium values in the pres- 
ence of strain, i.e., r) xyxy , f] xyzy and f]% zy are effective 
viscosities that have been renormalized by the relaxation 
of the director. 

Collecting, we obtain a set of equations for the dis- 
placement components that can be stated in the form of 
Eq. (|2.3[) with the components of the stress tensor given 
by 

Cxx(^) — C X xxx (^-0 ^xx H~ Cxxyyij^^} Uyy ~f~ ^xxzzi^') ^zz 
~\~ Cxxxz 

(id) U xz , (5.13a) 

(Jyy(uJ^J C X Xyy(j^} ^XX "I™ Cyyyy (Uj) ^yy H~~ Cyy ZZ [iU^ U ZZ 

+ C yyxz {uj) u xz , (5.13b) 

CZZ(^) = C XXZZ {Uj) U xx + Cyy ZZ [w) Uyy + C ZZZZ {uj) "zZ 

+ C zzxz (uj)u xz , (5.13c) 

<J V X{U) = \ C* yxy {uj) U X y + \ C^y yz {u) U yZ ~ B 1 dlUy 

- B 3 d x d 2 u y - 2 B 4 d 2 d z u y -2B 5 d 3 z u y , 

(5.13d) 

= \ C^ yyz {uj) u xy + \ C^ zyz {uj) u yz - B 2 dluy 

- B 3 d%u y ~2Bi dluy -2B 5 d x d z 2u y , 

(5.13e) 

(Txy(v) = \ C xyxy (u) U xy + \ C xyyz (uj) U yz 

-B 6 d^u x -B 8 d 3 y u z , (5.13f) 
cr 2y H = 2 C* yyz {uj) u xy + \ C* zyz {uj) u yz 

-B 7 d 3 u z -B s d 3 y u x , (5.13g) 

fflZ^) — O'zxi 1 -^) — 2 C XXXZ (U>) U xx + 2 C yyXZ {lO) Uyy 

+ \ C zzxz (u) u zz + \ C xzxz {lo) u xz , (5.13h) 

where Cijki (to) without superscript R is defined as in 
Eq. (|44|) and 

Cf x M = Q x -^r?f x . (5.14) 

Here, we have used a compact notation with indices £, \ 
running over the pairs xy and zy. We will on occasion re- 
turn to this notation below. When the shear moduli C^ x 
are nonzero and independent, in which case the bending 
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terms are unimportant for the hydrodynamic behavior, 
the stress tensor (|5.13l) is identical to the hydrodynamic 
stress tensor for an monoclinic solid. 



C. Sound velocities 

In order to determine the mode structure of SmC elas- 
tomers, we proceed as in our discussion of biaxial smec- 
tics by analyzing the sound-mode structure in the dis- 
sipationless limit. As before, any bending terms can be 
neglected in this analysis. Setting all viscosities to zero, 
we extract from the equations of motion the character- 
istic frequencies in the form of Eq. (|4.6|) . The result- 
ing sound velocities are plotted schematically in Fig. [5] 
There are 3 pairs of sound modes. One of these pairs 
(i) is associated with the soft deformations discussed in 
Sec. IV Al Its velocity vanishes for q along e y and e 2 so 
that when viewed in the z/z'-plane it has a clover leave 
like shape. The remaining 2 pairs are associated with 
non-soft deformations. In the incompressible limit, these 
pairs become purely transverse (ii) and longitudinal (iii), 
respectively. In the y'z' and a/z'-planes, their velocities 
are non-vanishing in all directions. Note that, since the 
velocity of pair (i) vanishes in directions where the veloc- 
ities of the other modes remain finite, SmC elastomers 
are, like nematic (54j and biaxial smectic leastomers, po- 
tential candidates for applications in acoustic polarizers. 




D. Mode structure in the incompressible limit 

Having found the general sound-mode structure in the 
nondissipative limit, we now turn to the full mode struc- 
ture in the incompressible limit. Due to the complex- 
ity of the equations of motion, we once again focus on 
the softness-related symmetry directions, cf. Fig. [5] (d) 
and (e). As in Sec. IIV C[ the equations of motion yield 
propagating modes, whose frequencies are of the form of 
Eq. (|4.8p and diffusive modes, whose frequencies are of 
the form of Eqs. (|4.9p or (|4.10[) . Specifics of the sound 
velocities, the diffusion constants and the bending terms 
entering these frequencies are given in the following. 

First, let us consider the case that q lies in the xz- 
plane. In this case the equation of motion for u y decou- 
ples from the equations of motion for u x and u z . The 
u y equation produces a pair of transverse propagating 
modes polarized along e y with a sound velocity 



C v = 1 r I cos9q x + sin6q z \ = \ — \q x >\ (5.15) 
4p y 4p 

and a diffusion constant 

Vxyxy Qx "f" ^ Vxyzy QxQ_z ~t~ Vzyzy Qz 



D 



p,y 



(5.16) 



FIG. 2: Schematic plots of sound velocities, (a), (b) and (c): 
spherical plots of respectively mode pairs (i), (ii) and (iii). (d) 
and (e): polar plots in the x'z'- and a;'z'-planes, respectively, 
of all 3 sound-mode pairs. Plot (a) is magnified relative to 
the remaining plots by a factor 2. 



In the soft direction, i.e. for q || e 2 <, the sound velocity C y 
vanishes and these propagating modes become diffusive 
with diffusion constant Dd.y = D p y and bending term 



By = 



B x qi 



B 



2Qz 



2B 4 q x q z + 2B 5 q x q 3 z 



(5.17) 



cf. Eq. (|4.9p . The equations of motion for u x and u z 
can be solved by decomposing (u x , u z ) into a longitudinal 
part ui along q and a transversal part ut along er = 
(q Zl 0, —q x ). In the incompressible limit ui vanishes. The 
equation of motion for ut produces a pair of propagating 
modes with polarization u || e^, with sound velocity 



C zzzz 2 C xxzz ]q x q z 



xxx ~r ^zzzz 

+ [C xxxz - C zzxz ]q x q z {ql ~ q 2 x ) 
+ \C xzxz {ql~ql)} 1/2 



(5.18) 
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and diffusion constant 



mulation with strain and director. 



?]x .1 



2^ 



2p 



xxzz -2 -2 



fxxxz 'fzzxz 



2p 



» - l -2 -2\ i Vxzxz i »2 ~2\2 



8p 



(5.19) 



Second and last, let us consider the case q |j e y . There 
is a pair of longitudinal propagating modes with u || 
that is suppressed in the incompressible limit. There is 
a pair of elastically soft diffusive modes polarized along 
e z i with 



A. Equations of motion 



B m , = 



ri R i 1 1 1 

Iz'y'z'y' 



sin 2 9 B 6 - sin 29 B 8 + cos 2 9 B 7 



(5.20a) 
(5.20b) 



where r$ , z , y , = sin 2 9 v^ yxy - sin 29 vg yzy + cos 2 9 vfc yzy . 
Finally, there is a pair of propagating modes polarized 
along e x < with 



C x i 



C/(A P ) . 
r, R , , , , 

fx y x y 



(5.21a) 
(5.21b) 



Now we will assemble effective equations of motion for 
the displacement variable u and the director. One of the 
main ingredient entering the general equations of mo- 
tion (|2.4j) is the elastic energy of the SmC phase. Here we 
need the elastic energy density with strain and director, 
Eq. (|5.ip . that comprises the stretching energy density 
/ u , Eq. (|5.3p . the coupling energy density / u , n > Eq. (|5.5[) . 
and the Frank energy density / n -Frankj Eq. (|5.7[) . The 
other important ingredient is the viscosity tensor. Again, 
we parametrize this tensor so that the entropy produc- 
tion density is of the form of Eq. (|5. 12|) , albeit now with 
rjijki replaced by Uijki (all unrenromalizcd thus far). 



with Vx'y'x'y' = cos 2 9 T/^ yxy + sin 29 r]* yzy 



Sin 2 0Vzyzy 



VI. SMECTIC-C ELASTOMERS - DYNAMICS 
WITH STRAIN AND DIRECTOR 

In this section we will treat the low-frequency, long- 
wavelength dynamics of SmC elastomers within our for- 



Recall that we have expanded the coupling and Frank 
energy densities about the equilibrium director (|5.2[) and 
that the deviation from the equilibrium director has only 
two components, which we denote for briefness by n x and 
n y . For convenience, we will in the following not work 
with n y directly but rather with the composite variable 
Q y defined in Eq. (|5.6[) instead. Then, the two equations 
of motion resulting from the director equation l|2.4a[) read 



[d t + TA + TK yy {W)]Q y +TK yx (W)n x = [A S d t - T A a]u xy + [X y/l - S 2 d t - T A 0\ u yz 

-TK yy (V){Sr] Ayx + V1-S 2 VAyz}, (6.1a) 
[d t + r K xx (V)}n x + r K xy {V) Q y = X S(l - S 2 ) d t [u xx - u zz ] + A (1 - 2S 2 ) v 7 ! - S 2 d t u xz 

-TK xy (V){S WAyx + a/ 1 ~~ S 2 TjAyz } + \/ 1 — S 2 dttjAxz j (6.1b) 



where 



^xs(^) ~~ ^xxxx $x R-xyxy $y ^xzxz $z ^ ^xxxz &x$z 3 (6.2a) 

^yy(^~) Kyxyx &x ^yyyy ^y -^-yzyz ^z ^ ^yxyz $x$z 1 (6.2b) 

■^xyi^J — [R-xxyy H~ -^xyyx) 9 x 0y [Kyy X z -^-xyyz) $y$z •> (6.2c) 

^yxi^J — [K-xxyy H~ ^xyyx) 9 x 0y [Kyy X z ~i~ -Kyxyz) $y$z • (6. 2d) 



The equations of motion resulting from the momentum density and displacement equations (|2.4b[) and (|2.4b[) can be 
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written in the form of Eq. (|2-3[) with a stress tensor given by 

fll(w) = C XXXX (uj) U xx + C XX yy(ljj) Uyy + C XXZZ (oj) U ZZ + C XXXZ (u>) U X 



+ A 5(1 - S 2 )K XX (V) n x + A 5(1 - S 2 )K xy (V) [Q v + Sr, Ayx + v 7 ! - S 2 VAyz] , 

&yyV*J) = C XX yy(uf) U XX -\~ Cyyyy{iS) Uyy ~t~ Cyy ZZ {iO^ U ZZ ~\~ Cyy XZ {ijj) U xz , 
CT ZZ (UJ) = C XXZZ (uj) U XX + Cyy ZZ {uj) Uyy + C ZZZZ {u) U zz + C ZZXZ {uj) U xz 

- A 5(1 - S 2 )K XX {V) n x - A 5(1 - S 2 )K xy (V) [Q w + S VAyx + y/l - 5 2 VAyz ] , 
CTyx(w) = 5 A [a + A 5]Q a + ^{a A [a + A 5] + C cos 2 9 + v xyxy d t }u xy 

+ |{/3 A [a + A 5] + Csin0 cos# + v yzyz d t }u yz , 
G yz {uj) = iA[/3 + AVl-5 2 jg y + i{aA[/3 + AVl-5 2 J + Csin<9 cos0 + v yzyz d t }u xy 

+ i{/3A[/3 + AVT^5 5 ] +Csm 2 6 + v yzyz d t }u yz: 
Cxy(u) = <7 yx (uj) + i (A - l)5iT 1/y (V){Qj / + 5 + V±--S 2 ~ VAyz} + 5 (A - l)5if H:c (V) 

Crrz^) = i Ccxx^Cj) flu + i Cy yxz (uj) U yy + „- C zza;z (cj) 1t zz + i C xzxz (uj) U xz 



+ [§ (A + 1) - A 5 2 ] \/l — 5 s " {tf ra (V) n x + ^(V) [Q„ + 5 ^ + VT^S 2 "^,] } , 

&Zx{t>}) = ^ C XXXZ (UJ) U xx + g Cyy XZ (iS) Uyy + ^ C ZZ£CZ (aj) U ZZ + ^ C XZ;( ; Z H IZ 

+ [i (A - 1) - A 5 2 ] a/I - 5^ {^(V) n x + K xy (V) [Q y + 5 r?^ + VT^VAyz] } , 
(J zy {oj) = cr zy {uj) , 



(6.3a) 
(6.3b) 

(6.3c) 

(6.3d) 

(6.3e) 
(6.3f) 

(6-3g) 

(6.3h) 
(6.3i) 



with Cijki (to) as defined in Eq. 



r 



Since these equations of motion are of considerable al- 
gebraic complexity, we will first consider the simplified 
case where all the K^i are set to zero. This will al- 
low us without too much effort to make contact to the 
hydrodynamic equations derived in Sec. IV Bl and to ex- 
tract renormalized elastic moduli that will be important 
for the behavior of SmC elastomers in rheology experi- 
ments. We will return to the full equations with the K^ki 
included further below. 

When the Kijki vanish, the dependence of the stress 
tensor and Eq. (|6.1ap on n x drops out. The latter equa- 
tion is then readily solved with the result 



Q y = -a 



1 + iu>T 3 



1 — iuiTi xy 1 — iuiT\ 
where we have introduced the relaxation times 

n = i/(TA), 



1 + iuJT 2 



(6.4) 



(6.5) 



r 2 = AVI- S 2 /{TA(3), t 3 = A5/(TAa). As we will 
see further below, our dynamical equations will predict 
non-hydrodynamic modes, which are essentially director 
modes, with a decay time ("mass") t\ implying that t\ 
is essentially the director relaxation time, t\ = r„. To 
estimate the value of r„, we note that the elastic con- 
stant A can be re-expressed in terms of the parameters of 
semi-microscopic models such as the neo-classical model 
of rubber elasticity [![. For details on the relation be- 
tween the elastic constants of our elastic energy (|5.ip 
and those of the underlying semi-microscopic models we 
refer to Ref. ■ Then, we use the known values of the 
semi-microscopic parameters to estimate the magnitude 



of A. We find A to be of the order of 10 3 Pa. We are 
not aware of any experimental or theoretical estimates 
of the rotational viscosity T _1 of SmC elastomers. We 
expect it to be somewhat larger than the values of about 
0.3 x 10 -1 Pas found in SmC liquid crystals [50f . perhaps 
r _1 lPas. This leads to r„ « 10~ 3 s as a rough esti- 
mate of the director relaxation time in SmC elastomers. 
This value is larger than the r„ ss 10 -2 s reported for ne- 
matic elastomers [HJ [H| , which is consistent with the ex- 
pectation that restoring forces on the director are greater 
in smectics than in nematics due to the layering. 

Feeding Eq. (|6.4[) into the stress tensor (|6.3p leads an 
effective stress tensor and corresponding effective equa- 
tions of motion in terms of the displacement variables 
only. This effective stress tensor has exactly the same 
form as the hydrodynamic stress tensor (|5.13p without 
the bending terms and with the renormalized frequency- 
dependent elastic moduli now given by 



ILUTl 



1 — IUT\ 

Q x _ ? ^f x + oV), 



■AA 



with renormalized viscosities 



'\x = v tx 



T A, 



£x 



(6.6) 



(6.7) 



In Eqs. (|6.6[) and (|6.7[) we have used our compact notation 
for index pairs established in Eq. (I5.14p . The individual 
A ix are given by A xyxy = a 2 {l+T 3 /Ti) 2 , A xyzy = a/3(l+ 
t 3 /t 1 )(1 + t 2 /t 1 ) and A zyzy = (3 2 (1 + t 2 /t 1 ) 2 . Note from 
Eq. (|6.6p that we can identify the bare viscosities rjijki 
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and Vijki as well as the renormalized viscosities r]^ and 

u^ x , indicating the consistency of our approaches without 
and with director. 



B. Mode structure in the incompressible limit 



Now turn to the equations of motion, Eqs. (|6.1[) 
and (|2.3|) in conjunction with (16. 3p . with the bend- 
ing terms included and considering the incompressible 
limit. As above, we will focus on the directions where 
Q = (fejO, q z ) and q = (0,^,0), respectively. Some of 
the steps involved in solving these equations of motion 
are sketched in appenix [B] . Contrary to the predictions 
of our purely hydrodynamical theories, the modes result- 
ing here are of three rather than of two different types. 
As was the case in the hydrodynamics of biaxial smec- 
tics and SmC's, there are propagating modes with fre- 
quencies in the form of Eq. (|4.8p and there are diffusive 
modes with frequencies as stated in Eqs. (|4.9p and (|4.10[) . 
In addition to these modes, however, our theory with 
strain and director produces non-hydrodynamic modes 
whose frequencies remain finite in the limit of vanishing 
wavevector. Their frequencies are of the form 



i D r . 



(6.8) 



These non-hydrodynamic modes with a zero-g decay time 
Ti are predominantly director relaxation modes, and T\ 
is the director relaxation time r n . In the following we 
will discuss the specifics of all three types of modes. 

Let us begin by considering case that q lies in the xz- 
plane. In this case the equation of motion for u y decou- 
ples from the equations of motion for u x and u z . The 
u y equation produces a set of transverse modes with u 
along e y . There is a non-hydrodynamic mode with 



D m ,y — 



Vxyxy Qx 



(6.9) 



Moreover, there are propagating modes with a sound 
velocity C y that is identical to the sound velocity C y 
resulting from our hydrodynamic theory of Sec. [Vj see 
Eq. (|5.15p . and with a diffusion constant 



v xyxy Qx ~l~ ^ v X yzy QxQz + V zyzy Qlz 

87 



(6.10) 



In the soft direction, i. e. for q |j e z r, these propagating 
modes become diffusive with Dd, y = D p y and 



B v = 



K x q% + K 2 qt + K 3 q'Ji + 2 K 4 q x q z + 2 K 5 q x q 



(6.11) 



where the K's are bending moduli that are combinations 
of the Frank elastic constants, the order parameter S as 



well as A, a, and (3. Our results for these K's are com- 
piled in appendix I A 31 We would like to emphasize, that 
the hydrodynamic theory and the theory with strain and 
director yield identical results for the sound velocity C y . 
Furthermore, Eqs. (|6.10p and (|5.16[) as well as Eqs. (|6.1ip 
and (|5.17p show that the results for the diffusion con- 
stants -Dp, j/ and the bending terms B y are in absolute 
agreement, when the identifications v xyxy = r] xyxy and 
so on as well as K\ = B\ and so on are made. 

The equations of motion for u x and u z are solved by 
decomposing (u Xl u z ) into a longitudinal part ui, that 
vanishes in the incompressible limit, and a transversal 
part ut- The equation of motion for ut yields a pair 
of propagating modes with a sound velocity as given in 
Eq. (|ET5)) . with 



xxxx H~ V Z zzz 2 Vxxzz ^9 ^2 

Txlz 



2p 



(6.12) 



v xxxz v zzxz . ~ / ~2 ~2\ , v -xzxz , ~2 -2\2 

+ o QxlAQz - 1x) + -5 — 9 2 - Qx) > 

2 p 8 p 

and with a u || er where er = {qz, 0, — q x ). Note the full 
agreement of Eqs. (|6.12p with Eq. (|5.19p when v^jki and 
Vijki are identified. 

Finally, we turn to q || e y . In this case here is a pair 
of longitudinal propagating modes with u |j e y that is 
suppressed in the incompressible limit. There is a non- 
hydrodynamic mode with 



Dn 



Vxyxy + v zyzy v zyzy 

4~P 



(6.13) 



where u lays in the rrz-plane with u x = cx(ti +t 3 )/[/3(ti + 
T2)] u z - There is a pair of diffusive modes with polariza- 
tion u || e z i that is related to elastically soft deforma- 
tions. The diffusion constant and the bending term of 
this mode pair reads 



D d}Z> = 



B z , = 



v H , , , , 

z'y'z'y' 

8p ' 

sin 2 9 K 6 - sin 29 K 8 + cos 2 9 K 7 



(6.14a) 



(6.14b) 



where vf; y , z , y , is the renormalized version of ^z'y'z'y', cf. 
Eq. (I6.19p . Last, there is a pair of propagating modes 
polarized along e x > with a sound velocity identical to that 
stated in Eq. (|5.21aP and with a diffusion constant 



D 



p,x' 



v R , , , , 

x' y' x'y' 



(6.15) 



cos 2 9 v* yxy - 



x29v R 



■sm 2 9v R yzy . 



wliere , , , 

x'y' x'y' ^ " xyxy 1 w xyzy 

Once more, we point out the full agreement of the dif- 
fusion constants and the bending terms obtained respec- 
tively through our purely hydrodynamical theory and our 
theory with strain and director when the proper identi- 
fications between the viscosities and between K e and B e 
etc. are made. This agreement on one hand signals the 
consistency of our two approaches. On the other hand, 
it reassures us that our algebra and our final results are 
correct. 
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C. Rheology 

In the hydrodynamic limit lot — ► 0, physical elastic 
moduli have the simple Maxwell (rubber-like) frequency 
dependence, i.e., C ijk i(u>) = C ljk i -iuv ijkl . At higher 
frequencies, they can exhibit more complex behavior aris- 
ing from non-hydrodynamic modes. Our simple theory 
focuses on one set of non-hydrodynamic modes, those 
associated with the director. It provides a good descrip- 
tion of the relaxation of the director, but it provides a 
good description of the dynamic elastic moduli only if r n 
is much greater than than any other relaxation time in 
the system, for example those associated with the relax- 
ation of polymer configurations. If we assume r„ ^> te, 
then our theory makes specific predictions about the fre- 
quency dependent elastic moduli and the corresponding 
experimentally relevant quatities, viz. the storage and 
loss moduli defined, respectively, as the real and imagi- 
nary parts of the frequency dependent moduli. 

One consequence of Eq. ()6.6p is is that the corre- 
sponding storage and loss moduli may exhibit uncon- 
ventional dips and plateaus when they are measured as 
functions of u>. Because rheology experiments typically 
probe a wide range of frequencies that reaches to fre- 
quencies well above the hydrodynamic regime, it seems 
useful to extend our discussion to frequencies exceeding 
lj <C Tg 1 . In oder to do so, we have to account for the fre- 
quency dependence of the viscosities and T. The Rouse 
model [53| , though it certainly does not provide a micro- 
scopically correct description of liquid-crystal elastomer 
dynamics, is known to provide for useful fitting curves 
for the storage and loss moduli of elastomers. This ob- 
servation leads us employ this model here and assume 
that the viscosities and T are respectively of the form 
Vijki(u) = v ljk if R (-iujT E ) and r(w) _1 = T' 1 f R (-iLOT E ) 
where 



3 

o 
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FIG. 3: Log- log plot of the reduced storage and loss moduli 
C'^ x (lj)/C^ x (solid lines) and C^ x {oj) /C^ x (dashed lines) ver- 
sus the reduced frequency lute as given respectively by the 
real and imaginary parts of Eq. (|6.17[) for t„/tb = 10 3 . For 
the purpose of illustration we have set, by and large arbitrar- 
ily- "{z/frsCfe) = 1 and AA (x /C (x = 3. 



experimentally. To discuss this implication, let us switch 
here to the rotated reference space coordinates x',y',z'. 
Experimentally, it will be hard to work with these coor- 
dinates, since the corresponding rotation angle 9 depends 
on temperature. Setting this problem aside, we note that 
there are two renormalized shear moduli that vanish in 
the limit to — > 0, viz. 



ium 

z'y'z' . t-±s±y'z>y'z' , 



1 — IL0T\ 



(6.18a) 



^ _ _ luJTl —AA 

L-\jTL X 'y'y> Z ' 3 



C x >yiy- Z >{u) = -lUVx'y'y'z' - — 



lLOT\ 



(6.18b) 



tt 2 x smh^^/x) 



{7T\/5cosh(7r-\/^) — sinh(7rv / ^)} , 

(6.16) 



is the well known Rouse function [53|. With these as- 
sumptions, we obtain the following phcnomcnological 
form for the complex moduli C^ x {lj) that incorporates 
the effects of both director and elastomer modes: 



a 



TE 



h R (-iu)T E ) 



(T n /TE) h R (-imT E ) 
1 + {Tn/TE) h R (-iLOTE) 



AA 



(6.17) 



where h R (x) = x f R (x). For r„ 3> te, the storage mod- 
uli C'^ x produced by Eq. (|6.17[) show an unconventional 
plateau and there is an unusual dip in the loss moduli 
C'J ' , cf. Fig. [3] If r„ » te, or r„ < te, this is not the 
case because elastomer modes will blur these plateaus 
and dips leading to conventional rubber-like behavior 
(not shown). 

There is another interesting theoretical implication of 
Eq. (|6.6p that, however, will be very difficult to detect 



with viscosities 



Vy' z'y'z' — sin Uxyxy sin 29 l^xyyz ~t~ COS 9 Vyzyz ■ 



(6.19a) 



Vx'y'y'z' = cos 29 v xyyz + | sin 29 (u 



'yzyz 1 xyxy I 



(6.19b) 



Ay'z'y'z' an d A x iyi y i z i can be inferred from Eqs. (|6.19p 
upon replacing v xyxy by A xyxy and so on. Since 

y'y'z'i^) van ish m the f° r lj ^ 0, 
we recover in this limit the ideal static soft elasticity dis- 
cussed in Sec. IV Al At non-vanishing frequency the sys- 
tem cannot be ideally soft but it can be nearly so for lj 
small. This deviation from ideal softness due to non-zero 
frequencies was first discussed in the context of nematic 
elastomers 5J], where this phenomenon was termed dy- 
namic soft elasticity. In a semi-soft SmC, the storage 
moduli C y , z , y , z ,(Lj) and G' x , yly , z ,{lj) in the rotated frame 
are nonzero at zero frequency. They will exhibit behavior 
similar to that Fig. [3] if r„ S> te ■ 
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VII. CONCLUDING REMARKS 

In summary, we have studied the low-frequency, long- 
wavelength dynamics of Srm4, biaxial smectic and SmC 
elastomers, assuming that these materials have been 
crosslinked in the SmA phase. We employed two different 
but related approaches: one formulation that does not 
explicitly involve the Frank director and that describes 
pure hydrodynamics and a second formulation that fea- 
tures the director and that describes slow modes beyond 
the hydrodynamic limit. 

The hydrodynamics of SmA elastomers is qualitatively 
the same as that of conventional uniaxial rubbers. Be- 
yond the hydrodynamic regime, however, the director has 
an impact on the dynamics of SmA elastomers, in which 
the genuine difference between SmA and conventional 
uniaxial elastomers shows up. The low-frequency, long- 
wavelength dynamics of SmA elastomers, as described by 
our theory with the director, is qualitatively the same as 
that of nematic elastomers, up to those aspects, where 
the value of the modulus Ctf for shears in the planes 
containing the director enters. For example, nematic 
and SmA elastomers possess 3 pairs of hydrodynamic 
sound modes. One of the qualitative differences is, that 
a soft nematic elastomer, where Cf = 0, has two pairs of 
sound modes whose sound velocity vanishes in the soft- 
ness related symmetry directions, whereas for a SmA 
elastomers, where C§ is significantly larger than zero, 
sound velocities are non-vanishing in all directions. The 
other qualitative difference is that nematic elastomers 
exhibit dynamic (semi-)soft elasticity, i.e., the general- 
ization of (semi-)soft elasticity to non-zero frequencies. 
SmA elastomers crosslinked in the SmA phase, since they 
are not soft or semi-soft, cannot have this property. 

Biaxial smectic and SmC elastomers possess, like ne- 
matic and SmA elastomers, three pairs of hydrodynamic 
sound modes. One of the sound- mode pairs in either elas- 
tomer is associated with soft deformations. The sound 
velocity of these pairs, which become purely transverse in 
the incompressible limit, vanishes in the softness-related 
symmetry directions. The remaining two sound sound 
mode pairs are associated with non-soft deformations. 
Their velocities remain finite in the directions in which 
the velocities of the softness-related modes vanishes. In 
the incompressible limit, these mode pairs become purely 
transverse and longitudinal, and hence the latter are ef- 
fectively suppressed. The vanishing of the velocities of 
the softness-related modes in certain directions could be 
exploited, in principle, to separate these modes from the 
remaining modes, which have different transverse polar- 
izations than the softness-related modes. Therefore, at 
least theoretically, soft biaxial smectic and SmC elas- 
tomers could be used to polarize acoustic waves. 

Our theory with strain and director reveals that SmC 
elastomers may have an unconventional behavior in rhe- 
ology experiments. It predicts that, as in nematic and 
SmA elastomers, certain storage and loss moduli may 
show, respectively, a plateau and an associated dip in 



the frequency range r" 1 < w < t^ 1 . It also predicts 
that soft samples of these materials exhibit dynamic soft 
elasticity. Recently there has been a controversy whether 
this unconventional rheology and dynamic soft elasticity 
can be observed in liquid crystal elastomers or not. Thus 
far, this controversy revolved around nematic elastomers. 
The few samples of nematic elastomers on which rheol- 
ogy measurements have been performed, produced lit- 
tle evidence for unconventional rheology or dynamic soft 
elasticity. However, due to the limited amount of data 
that is currently available and due to differences in their 
interpretation, the experimental picture still gives reason 
for debate; both phenomena are neither clearly verified 
nor clearly ruled out. The idea to use SmC elastomers as 
an alternative testing ground for these phenomena seems 
appealing. As far as dynamic soft or semi-soft elasticity 
is concerend, however, it will be difficult if not impossible 
to realize soft or semi-soft oszillatory shears experimen- 
tally, because one has to use a very specific, temperature- 
dependent coordinate system. Regarding the plateau 
and the dip in, respectively, the storage and loss mod- 
uli, these are less likely to exist physically in smectics 
than in nematics, because the director relaxation time in 
smectic elastomers is probably shorter than in nematic 
elastomers. 

Despite these difficulties, we hope that our theory en- 
courages rheology experiments on smectic elastomers, in 
particular, to see if unconventional rheology exist or not. 
Moreover, we hope that our work motivates experimental 
investigations of sound velocities in smectic elastomers, 
for example Brillouin scattering experiments. 
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APPENDIX A: ELASTIC MODULI 



This Appendix collects specifics about the Frank elas- 
tic constants defined by Eq. (|5.7p . the bending moduli 
defined by Eq. (|5.10|) and the effective bending moduli 
appearing in Sec. IVI Bl 
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1. Frank elastic constants appearing in Eq. (|5.7p 
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APPENDIX B: SOLVING THE EQUATIONS OF 
MOTION IN THE SMC* DYNAMICS WITH 
STRAIN AND DIRECTOR 



Here we sketch our analysis of the equations of motion, 
Eqs. (|6.1|) and (|2.3p together with (|6.3|) . in order to obtain 
the modes discussed in Sec. IVIBI 



Let us start be introducing the Fourier transforms of the differential operators K XX (V) and so on defined in Eq. (|6.2 



q K xx = K, 



v.i: Qx ^ xyxy Qy 



K xz 



2K 



xxx z QxQz ) 



Q Kyy — K yxyX q x + Kyyyy Q y + K yzyz + 2 K yxyZ q X Qz , 

q 2 K xy 



Qy{(.^-xxyy ~\~ -^-xy 



ryyx) Qx 

Q = Qy{[\K XX yy -\- K xyyx) Qx 



(,-^yyxz -K-xyyz) Qz\ 

{-KyyxZ H~ Kyxyz} Qz^ 



(Bla) 
(Bib) 
(Blc) 
(Bid) 



When q y — or q x 



0, the cross terms q 2 K xy and q 2 K yx vanish. Writing Eqs. (|6.ip in momentum space 



18 



one observes, that in this case, the equations for Q y and n x , decouple. This decoupling simplifies the solution of the 
equations of motion considerably and we will limit our following consideration to momenta where this simplification 
applies. 

When Eqs. (|6.1[) decouple, they are readily solved for given displacements u x , u y and u z with the result 
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Substituting Eqs. (|B2| into the stress tensor (|6.3| . one obtains effective equations of motion for the displacements 
only. 

As a specific example, let us now consider the case q y = in some detail. The second case, q x = q z = 0, can be 
treated by similar means and will be left as an exercise to the reader. For q y = 0, the effective equation of motion for 
u y decouples from those for u x and u z , which remain coupled. Expanding in powers of q and oj, we find the equation 
of motion for u y to read 

P Uy = | | C^yiUj) q^l C^yjuj) q^ + j CKyy^Uj) ^ 

+ Kx qi + K 2 q\ + K 3 q 2 x q 2 z + 2 K± q 3 x q z + 2 K 5 q x q 3 2 } U y , (B3) 

where we have dropped higher order terms that do not affect our results. The bending moduli K\ and so on are 
as defined in App. IA 31 Solving Eq. (|B3[) for frequencies in the form of a power series in q then leads to the non- 
hydrodynamic mode with the diffusion constant stated in Eq. (|6.9|) and the hydrodynamic modes with sound velocities, 
diffusion constants and bending contributions as given in Eqs. (|5.15p . (|6.10[) and (|6.1ip . 

Finally, let us look at the coupled equations for u x and u z . Proceeding like above, we obtain 
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as the leading contributions. Equations (|B4[) can be decoupled in the incompressible limit. To this end, we switch now 
from the coordinate system with basis {e x ,e yi e z } to a rotated system with basis {q,e y ,eT} with = (q z , 0, — q x ). 
In the incompressible limit the longitudinal component U[ of the displacement vanishes. The equation of motion for 
the transversal component ut along follows from Eqs. (|B4[) as 



(B5) 
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Solving this equation for lo in the form of a power series in q results in the propagating modes with a sound velocity 

and a diffusion constant as resented in Eqs. (|5. 18|) and (|6.12|) . 
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